Mechanism and Lineshapes of Sub-Fourier Resonances 
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Subtle internal interference effects allow quantum-chaotic systems to display "sub-Fourier" reso- 
nances, i.e. to distinguish two neighboring driving frequencies in a time shorter than the inverse of 
the difference of the two frequencies. We report experiments on the atomic version of the kicked 
rotor showing the unusual properties of the sub-Fourier resonances, and develop a theoretical ap- 
proach (based on the Floquet theorem) explaining these properties, and correctly predicting the 
widths and lineshapes. 
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A major difference between classical and quantum sys- 
tems is the existence in the latter of interferences between 
various paths. At long times, there are typically a large 
number of trajectories and one could expect that, in the 
average, the contributions of the various paths have un- 
correlated phases, leading to vanishingly small interfer- 
ence terms; in this view, quantum and classical transport 
should be identical. This simple expectation is however 
too naive. Even for classically chaotic systems where tra- 
jectories are very complicated and proliferate rapidly, it 
is not true that the phases of the various contributions 
are uncorrelated. A dramatic example is the so-called dy- 
namical localization (DL) taking place in a periodically 
kicked rotor, described by the following Hamiltonian: 
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where 9 is the 27r— periodic position of the rotor, p the 
conjugate momentum, T the period and K is propor- 
tional to the strength of the kicks). For if > 5, the 
classical dynamics is essentially chaotic. Although deter- 
ministic, the motion can be, on the average, described 
as a diffusion in momentum space, with the average ki- 
netic energy growing linearly with time. The situation is 
completely different for the quantum evolution: after an 
initial diffusive growth, the average kinetic energy satu- 
rates after some break-time. 

DL, leading to a complete freeze of the classical diffu- 
sion, is a non-trivial effect, as it means that the quantum 
interferences between classical diffusive trajectories are 
on the average completely destructive at long time. A 
fruitful insight is obtained by using the Floquet theorem. 
The so-called Floquet states \ipk) are defined as the eigen- 
states of the unitary evolution operator U{T) associated 
with Hq over one period: U{T)\ipk) = exp (— iefc) |<pfc) 
where is the eigenphase. The temporal evolution of 
any state after n periods is: 

l^(nT)) = [U{T)r\m) =Y.<'k 6-^"^'= Wk) (2) 



with Cfe = 1-0(0)). The evolution of any quantity can 
be calculated using the Floquet basis, for example: 

{p^{nT)) =J2ckcl, e-»(^'=-^''') iipk'\p^\iPk). (3) 

k,k' 

The Floquet states \ipk) of the chaotic kicked rotor are 
well known: they are on the average exponentially lo- 
calized in momentum space around a most probable mo- 
mentum pk , with a characteristic localization length £ Q . 
Such a localization - at the origin of DL - is far from 
obvious and is closely related to the Anderson localiza- 
tion in time-independent disordered one-dimensional sys- 
tems i2j- As it is typical for classically chaotic systems, 
the Floquet spectrum is highly sensitive to changes of 
the parameters, and the momentum wavefunction ipk(j>) 
displays large fluctuations. 

If the initial state |V'(0)) is well localized in momentum 
space around, say, zero momentum, Eq. ^ implies that 
only Floquet states with roughly \pk\ ^ i ("important" 
Floquet states) will play a significant role in the dynam- 
ics. Eq. (jSJ is a coherent sum over Floquet states. How- 
ever, as times goes on, non-diagonal interference terms 
accumulate larger and larger phases. In a typical chaotic 
system, these phases will be uncorrelated at long times, 
and interference terms will on the average cancel out, 
leading to an incoherent sum: 
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This equation is valid when DL is established. How long 
does it take for the phases n{ek — ek' ) to be of the order of 
27r? This can be simply estimated from the level spacing 
between important Floquet states and turns out to be 
roughly tbreak = £T , while (p^) saturates to a value oc P. 

DL has been observed in experiments with cold atoms 
exposed to kicks of a time-periodic, far detuned, standing 
laser wave H, ^ |B| - Our experimental setup is described 
in 0]. Basically, cold cesium atoms are produced in a 
standard magneto-optical trap. The trap is turned off, 
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FIG. 1: Sub-Fourier resonances. The measured atomic pop- 
ulation in the zero-momentum class, after 20 (upper solid 
curve) and 100 (lower solid curve) double kicks is plotted 
vs. the ratio r of the periods. The period of the first sequence 
is T = 27.8/is, and /C' ~ 10. At r = 1, dynamical localization 
is responsible for the large number of zero-momentum atoms. 
Away from the exact resonance, localization is progressively 
destroyed. Note the narrowness (5 times smaller than Fourier 
limit) and the triangular shape of the resonance line. Dashed 
lines are fits to the experimental lines using Eq. JHJ. 
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and pulses of a far-detuned (9.2 GHz ^ 1700 F) standing 
wave (around 90 mW in each direction) are applied. At 
the end of the pulse series, counter-propagating phase- 
coherent Raman beams perform velocity-selective Ra- 
man stimulated transitions between the hyperfine ground 
state sublevels Fg = 4 and Fg = 3. A resonant probe 
beam is used to estimate the fraction of transfered atoms, 
thus measuring the population of a velocity class. Re- 
peated measurements allow to reconstruct the atomic 
momentum distribution P{p). 

When the system is exposed to a two-frequency non- 
periodic driving, no DL is expected 0- Consider the 
Hamiltonian: 
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e ^ [5{t - nT) + 5{t - nrT - XT)] . 
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where r is the ratio of periods and A the initial phase be- 
tween the two kick sequences. If r is rational, the system 
is strictly time periodic and DL takes place, but DL is 
rapidly destroyed around any rational number, as exper- 
imentally observed in 0. Fig. H] shows the experimen- 
tally measured population in the zero-momentum class 
P{p = 0) after 20 and 100 double kicks as a function 
r. It displays a spectacular peak at r = 1 and a sharp 
decrease on both sides, due to the destruction of DL. 
There are two surprising features: (i) The resonance is 
very narrow: after N kicks, it could be argued that the 
two quasi-periods can be distinguished only if they differ 
by 1/N (in relative value). This would predict a width of 
the order of Ar — 1/N — 0.01 for 100 kicks, whereas we 
experimentally observe 0.0018 {ii) This "sub- Fourier 
resonance is not smooth, but has a marked cusp at the 



FIG. 2: (a) Averaged squared momentum {p^) {p is in recoil 
momentum units) vs. number of double kicks. All curves show 
DL after 5-10 kicks. After this break-time, the dynamics is 
frozen for the periodic system (r=l, solid curve), but a resid- 
ual quantum diffusion is observed in the quasi-periodic case, 
for r = 1 + 0.00111 and r = 1 - 0.00111 (dotted curves) and 
r — 1 ii 0.00222 (dashed curves), (b) The quantum diffusion 
constant (in units of the classical diffusion constant) versus 
r, from data shown in (a). It (almost) vanishes at r = 1 
because of DL, and rapidly increases symmetrically on both 
sides, with a triangular cusp. 



maximum. The aim of this paper is to discuss the physi- 
cal mechanism responsible for the destruction of DL and 
to give an explanation for these two unexpected features. 

Let us now consider the system evolution. Fig. ^a,) 
shows (p^) as a function of time, for various r values. 
Note that we measure P{p ~ 0) and then deduce (p^) as- 
suming that (p^) oc l/[P{p = 0)]^. For r = 1, DL is visi- 
ble after a break-time of the order of 5-10 periods, as (p^) 
saturates to a finite value. For the other values of r, the 
short-time dynamics is similar, but, for t > tbrcak, ip'^) in- 
creases slowly, linearly with time. In other words, in the 
vicinity of the resonance, the quantum transport is not 
completely frozen, a residual quantum diffusive process 
takes place, with a diffusion constant much smaller than 
the classical constant observed at short times. Fig. |2Jb) 
shows the quantum diffusion constant, as a function of 
r. It displays a sharp minimum at r = 1 with a char- 
acteristic triangular shape. Due to residual spontaneous 
emission by the kicked atoms ,9], the quantum diffusion 
constant at r = 1 is not strictly zero. 

The key point to understand the dynamics of the sys- 
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FIG. 3: Eigenphases of the evolution operator over one 
period for the doubly-kicked rotor, corresponding to the ex- 
perimental conditions, versus the phase A. The finite dura- 
tion of the kicks (800 ns) is taken into account. We have 
only plotted Floquet states having a significative weight: 
\{ip{0)\(pk{)^o))f > 10"" (thin line), or > 10"^ (thick line). 
States appear (disappear) as their weights go above (below) 
the threshold. Note that Floquet states rapidly change when 
avoided crossings with other states are encountered. 



tern is to realize that, near r = 1, the sequence of kicks is 
very similar to a periodic doubly- kicked rotor for r = 1, 
except that the phase A between the two sequences slowly 
drifts. The evolution operator U of the quasi-periodic 
system from time nT to {n + 1)T is thus given by the 
evolution operator U of the periodic doubly-kicked ro- 
tor for slowly changing A: U{r; Xo;nT ^ {n + 1)T) = 
U{Xq + n{r — 1)), where Aq denotes the initial phase be- 
tween the two pulse sequences. The total evolution op- 
erator from t = to {N — 1)T can thus be written as the 
product rin^o^ C^(^o + 'n{r — 1)). For small enough |r— 1|, 
the adiabatic approximation can be used, which states 
that, if the system is initially prepared in a Floquet eigen- 
state of U{Xo) it remains in an eigenstate of the instan- 
taneous evolution operator U{X), following the continu- 
ous deformation of the Floquet eigenstate of the periodic 
doubly-kicked rotor. This is illustrated in Fig. which 
shows the Floquet spectrum of the periodic doubly-kicked 
rotor, obtained by numerically diagonalizing U{X). When 
A is varied, the energy levels evolve in a complex manner. 
The complicated "spaghetti" behaviour, characteristic of 
quantum-chaotic systems, with a large number of avoided 
crossings (AGs) is directly responsible for the sub-Fourier 
character of the resonance lines, as we show below. 

In our case, the initial state is a linear combination of 
Floquet eigenstates. In the adiabatic approximation, it 
remains a linear combination of the instantaneous Flo- 
quet states with the same weights (the phases change, 
but the squared moduli remain constant). As discussed 
above, the coherences between Floquet states play a role 
only at short time, and interference terms average to zero 



after the break-time, Eq. This implies that Eq. @ 
remains valid close to the resonance, provided one uses 
the instantaneous Floquet eigenbasis: 

{p'w) « E \{m\vki\o))\' {M^)\p'\v'k{x)). (6) 

k 

The weights |(?/'(0)|(^fe(Ao))P are estimated at the initial 
time (A — Aq,), while the average values {ipk{X)\p'^\(pki^)) 
are calculated at the final value A = Ao-l-iV(r— 1). Eq. © 
explains the experimental observations. The fundamen- 
tal point is that Eq. © describes a correlation between 
the localization properties of the Floquet eigenstates for 
two values of A. The \{ip{0)\'Pk{Xo))\'^ term gives impor- 
tant weights to the Floquet states localized near zero mo- 
mentum. As A moves away from Aq, the important Flo- 
quet states evolve and, on the average, their center moves 
away from zero momentum. Thus, one expects (p^(A)) 
to have a minimum at A = Aq and to rapidly increase on 
both sides of Aq. Since A — Aq is proportional to r — 1, (p^) 
shall present a sharp minimum at r = 1; simultaneously, 
the population in the zero-momentum class shall present 
a sharp maximum at r = 1, which is experimentally ob- 
served in Fig. n 

What determines the linewidth of the resonances? A 
state \(pk(X)) evolves because of its AGs with other Flo- 
quet eigenstates. Tiny AGs may be crossed diabatically, 
leading to a large modification of the state. On the av- 
erage, the scale AA = A — Aq on which \(pk{X)) loses the 
localization property of \ifk{Xo)) is the distance to the 
next AG. We immediately deduce that the width A/ of 
the sub- Fourier line is such that: 

NTAf = 2AA. (7) 

Eq. lO also explains why chaos is necessary to ob- 
serve sub-Fourier resonances. If the classical dynamics 
is regular, the Floquet eigenstates evolve smoothly with 
the parameter A; a change in A of the order of one is 
thus required to significantly modify the Floquet states: 
2AA ~ 1, which is the Fourier limit. On the contrary, in a 
classically chaotic system, the level dynamics is known to 
be complex with plenty of AGs (Fig.lSJ, A A is then much 
smaller than unity, leading to highly sub-Fourier reso- 
nances (up to a factor 37 observed in j^])- Eq. also 
predicts the linewidth to be inversely proportional to the 
temporal length of the kick sequence beyond the break- 
time (i.e. the sub- Fourier character is independent of iV), 
as numericaly observed 0. The value of A A depends on 
the detailed dynamics of the system. It can be roughly 
estimated by visual inspection of the quasi-energy level 
dynamics. Fig. |21 The level spacing between important 
Floquet states is the inverse of the number of such states, 
itself roughly equal to the localization length, of the order 
of 5 for the present set of parameters. The typical rela- 
tive slope of neighbor levels in Fig. |3|is of the order of 4, 
thus predicting AA of the order of 0.05 and a "factor 10" 
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sub-Fourier line, about twice the experimentally observed 
factor (the difference is due to the transverse profile of 
the laser mode leading to spatial inhomogeneities in K). 

From Eq. Q and the previous analysis, it is expected 
that (p^) significantly increases from its minimum value 
at A = Ao as A varies to Aq + AA. It seems reasonable to 
assume that this increase is initially quadratic in A — Aq 
[l(ij | . The situation here is more complicated because Flo- 
quet states are localized in momentum space. AGs be- 
tween Floquet states localized around the same momen- 
tum, will be rather large (this is what determines AA) 
whereas AGs between states localized a distance L ^ i 
apart in momentum space are typically much smaller and 
scale like exp{—L/£). There is thus a very broad distri- 
bution of AG widths, with a large number of tiny AGs. 
A tiny AG typically extends over a small A interval and 
thus tends to produce small values of AA. The increase 
of (p^) with A thus depends on the number of small AGs 
encountered. In the presence of exponential localization, 
the AG density scales with size C as 1/C for C — > 0, and 
(p^(A))-(p^(Ao)) shall behave like |A-Ao| [ui, producing 
the cusp experimentally observed in the resonance line. 
Fig. ^ and in the diffusion constant. Fig. |2 The large 
number of extremely small C is responsible for the singu- 
larity of the sub- Fourier resonance line. Another conse- 
quence is the diffusive behaviour observed in the vicinity 
of the resonance, see Fig. |2b). Indeed, (p^(A)) — (p^(Ao)) 
increases linearly with | A — Aq | , itself a linear function of 
time and |r — 1|. Thus, our model correctly predicts two 
non-trivial properties : (p^) increases linearly with time 
and the corresponding diffusion constant is proportional 
to \r - 1|. 

There will always be some degree of nonadiabaticity. 
Whatever small |r— 1 1 is, tiny enough AGs will be crossed 
diabatically. This puts a lower bound on the size of the 
AGs effectively participating in the quantum dynamics 
of the system and produces a rounding of the top of the 
sub-Fourier resonance line, which is too small to be seen 
in the experiment after 100 kicks, but is easily visible 
after 20 kicks, Fig. [2 

Our approach concentrates on the immediate vicinity 
of the resonance. What happens in the wings of the sub- 
Fourier line? Eq. indicates that this depends on the 
residual correlation between |(/3fe(Ao)) and \ipk{X)) for |A — 
AqI > AA. A quantitative answer to this question cannot 
be given, as there is no typical AG width. However, it 
seems clear that the quantum diffusion constant does not 
exceed the classical one, which correspons to the killing 
of all interference terms. Random Matrix Theory tells us 
that this type of parametric correlation usually decays 
algebraically with A — Aq. We thus propose the following 
anzatz: 



the saturation value of due to DL. This equation fits 
(using P{p = 0) cx {p'^)~^^'^) very well the experimen- 
tal curves in Fig. ^ reproducing the linear behaviour at 
the center of the resonance and the classical diffusion in 
the wings. Nevertheless, because of the laser inhomo- 
geneities, the parameters are fitted, not extracted from 
the previous analysis. 

In summary, we have developed a theoretical approach 
for the mechanism of sub- Fourier resonances, which cor- 
rectly predicts the unexpected observed features. Devia- 
tions from exact periodicity are treated in the framework 
of the adiabatic approximation for the Floquet spectrum 
of the system, which thus goes beyond a pertubative ap- 
proach. The dynamics is governed by instantaneous Flo- 
quet eigenstates, which are non-trivial objects, as they 
need to have very well defined internal (i.e. between the 
various parts of the wavefunction for each Floquet state) 
phase coherence to be stationary states of the strictly 
periodic system, but at the same time, the inter-state 
coherences do not play any role (beyond the break-time, 
Floquet states effectively add incoherently). We are not 
aware of any other quantum system where the dynam- 
ics is dominated by such an incoherent sum of internally 
extremely coherent states. It shows that the role of inter- 
ferences in quantum mechanics is far from obvious, and 
can produce unexpected behaviors. 
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nologies de Lille. Laboratoire Kastler Brossel is labora- 
toire de I'Universite Pierre et Marie Gurie et de I'Ecole 
Normale Superieure, UMR 8552 du GNRS. GPU time on 
various computers has been provided by IDRIS. 
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|r- 1| -I- AA/tbroak 

where Dd is the classical diffusion constant and {p'^)dl 
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